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Abstract 

We compute the one-loop dispersion relations at finite temperature for 
quarks, charged leptons and neutrinos in the Minimal Standard Model. 
The dispersion relations are calculated in two different plasma situations: 
for a vacuum expectation value v of the Higgs field u / (broken elec- 
troweak symmetry) and for v = (unbroken electroweak symmetry). The 
flavour and chiral non-degeneracy of the quasi-particle spectrum is studied. 
Numerical results show that the thermal effective masses for fermions in the 
broken phase have a smaller value than those in the unbroken phase. The 
temperature dependence of the top quark and electron neutrino thermal 
effective masses is also presented. Gauge invariance of one- loop dispersion 
relations is studied. 
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1 Introduction 



The behaviour of the plasma in the high temperature regime plays a fundamen- 
tal role in the attempts to explain several puzzles in cosmology, as for example 
the present baryon asymmetry of the universe [|l], in the electroweak scenario. 
When the temperature of the universe (T) was near the critical one of the cos- 
mological electroweak phase transition (Tc), the dynamics of the plasma was 
governed by the interactions of the Minimal Standard Model (MSM). It is now 
well known that the interaction of a fermion with the thermal background in a 
plasma modifies the poles of the fermion propagator with respect to those at tree 
level. In a pioneering work ||, Weldon studied the modification of the fermion 
dispersion relation, i.e. poles of the fermion propagator, due to high temperature 
effects for a chirally invariant gauge theory (for instance, QCD or QED with 
massless fermions). In the same work, he showed that for a chirally invariant 
gauge theory with parity violation, as is the case for the SU(2)l <S> U(1)y model, 
the dispersion relations for the right- and left-handed massless fermions are in- 
dependent. This decoupling implies that the thermal effective fermion masses, 
defined as the energy values for zero momentum, are different. 

For a parity preserving gauge theory at finite temperature, the fermion disper- 
sion relation has two possible solutions. These solutions (branches in the phase 
space), both of positive energy, describe the propagation of the fermionic excita- 
tions of the plasma (quasi-particles) through the thermal medium. In the litera- 
ture, these solutions are generally known as normal and abnormal branches. Each 
point of a branch corresponds to the energy and momentum values of a quasi- 
particle. The quasi-particles described by the normal branch (quasi-fermions) 
have the helicity equal to their chirality, while those described by the abnormal 
branch (quasi-holes) have the helicity opposite to their chirality. This was stud- 
ied in detail by Weldon for massless quarks in QCD {|, || and by Pisarski for 
light massive fermions in QCD/QED |J. The normal branch is so called be- 
cause it becomes the ordinary fermionic branch at large momentum (large mass), 
k (Mf) > gT, where g stands for the coupling constant and Mf for the fermion 
mass. The abnormal branch owes its name to the fact that it appears as an 
additional physical solution only if temperature effects are taken into account. 
Quasi-particles reflect a collective behaviour of the plasma system at low mo- 
menta 0. The collective behaviour of the quasi-holes only exists for k < gT, 
where their creation probabilities are non-zero. Moreover, the QCD/QED quasi- 
holes are physically possible for Mf < gT and they disappear from the spectrum 
forM/~T[§|. 

The theoretical expectation that a cosmological electroweak phase transition 
H took place during the cooling of the early universe is motivated by the existence 
of the electroweak spontaneous symmetry breaking (ESSB) in the MSM and its 
extensions. The idea of the electroweak phase transition in the framework of 
the MSM allows the consideration of the two following phases. For T > Tc, the 
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gauge symmetry of the fermion fields is SU(3)c® SU(2)l<S>U(1)y, v = and the 
plasma scenario is the electroweak unbroken phase. For T < To, the electroweak 
symmetry SU(2)l <g> U(l)y is spontaneously broken. The gauge symmetry of 
the fermion fields is SU(3)c <8> U(1)em, v ^ 0, and the plasma scenario is the 
electroweak broken phase. 

Fermionic excitations are present in both electroweak phases of the plasma. 
The spectrum of quasi-particles can be obtained from the thermal loop contri- 
bution to the fermion self-energy. The energy of the quasi-particles is of order 
g s T. Quasi-particles, during their propagation in the plasma, scatter with the 



thermal medium and this phenomenon induces a damping rate 7 ||10|| . This 7 
is proportional to the imaginary part of the self-energy |TT[| . The energy of the 
quasi-particles compared to 7 is large and it is thus possible to speak about the 
existence of quasi-particles as coherent excited states. Quasi-particles have a fi- 
nite lifetime ~ I/27 and they can eventually turn into a new state, which is out 
of phase with the initial one. If the group velocity is \v\, the mean free path is 
\v\/2 r y. The quantum spatial coherence, i.e. the phase relation between points 
separated by a distance > \v\/2^, is lost p], [12]. Since we are interested in ob- 
taining the branches which describe the quasi-particles, only the real part of the 
one-loop contribution to the self-energy will be calculated. 

The main purpose of the present work is to compute exactly the one-loop 
dispersion relations for the different fermion sectors of the MSM for the following 
situations of the plasma: for the electroweak broken phase and for the electroweak 
unbroken one. We study numerically the behaviour of the normal and abnormal 
branches taking into account all terms in the fermion self-energy^. The flavour 
and chiral non-degeneracy of the quasi-particle spectrum is studied. Our atten- 
tion is focused in the low-momentum regime where collective phenomena develop. 
Since mass and T corrections are considered numerically, we obtain a novel re- 
sult, namely that thermal effective masses for fermions have a lower value in the 
broken phase than in the unbroken one. The existence of a similar shift was 



observed in |TJ] where, through a lattice calculation, it was found that thermal 
effective masses for massive bosons are smaller than those for massless bosons. 
We also confirm the particular behaviour of the branches in the broken phase 
found in |J : at variance with what happens in the unbroken phase, an energy 
gap of order Mf is obtained between the two branches that now have hyperbolic 
form. This behaviour is crucial for the mechanism of electroweak baryogenesis 
produced at the scattering of quasi-particles off the boundary created during the 
first order electroweak phase transition. 

The T dependence for the top quark dispersion relations in the broken phase 
and for the neutrino dispersion relations in both phases is also studied. We 



3 By exact numerical computation of the fermion dispersion relations, we mean that we 
have included the leading and non-leading T corrections and the fermionic and bosonic mass 
corrections. 
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have found, for the top quark, that the thermal effective mass for two of the 
four branches can have negative value for T < M t . In order to study the gauge 
independence of the dispersion relations, we compare one-loop computations in 
Feynman and Landau gauges. The dispersion relations are only gauge indepen- 
dent at leading order. 

In section 2, we define the fermion dispersion relations for the different types of 
gauge theories involved in the calculations. In section 3, the generic gauge boson 
and scalar diagrams of Fig. 2 are calculated in the Feynman gauge and the real 
part of the finite temperature contribution to the fermion self-energy is given. 
In section 4, the Lorentz-invariant functions for quarks, charged leptons and 
neutrinos in the MSM are calculated and the corresponding dispersion relations 
in both phases are derived. In section 5, we study in detail the gauge invariance 
for the one-loop dispersion relations for massless and massive quarks in QCD. 
The results are extended to the MSM dispersion relations. Numerical results of 
the dispersion relation computations are presented in section 6. Our conclusions 
are summarized in section 7. Results about the calculation of the real part of 
the fermion self-energy in the Landau gauge are presented in Appendix B. In 
Appendix C, the Lorentz-invariant functions calculated in the Landau gauge are 
also presented. 

2 Some useful definitions 

The fermionic dispersion relations, for the different cases in which the chiral 
and parity symmetry are preserved or violated by a gauge theory, are presented 
in this section. The fermionic spectrum of the quasi-particles in a plasma is 
described by the branches obtained from the fermion dispersion relation, which 
is derived from the real part of the thermal contribution to the self-energy. The 
fermion self-energy is written as a linear combination of the matrices 1, Ufa t/i and 
where is the four- velocity of the plasma, with u^u^ = 1, and is the 
fermion momentum. In the one-loop calculation, Ififi terms are not generated and 
the real part of the thermal contribution to the fermion self-energy, for a gauge 
theory with parity and chirality violations, has the following form [[|: 

Re E'{K) = -${a L L + a R R) - i(b L L + b R R) - Mj(c L L + c R R), (2. 1) 

where L = |(1 — 75) and R = |(1 + 75) are the left- and right-handed chiral 
projectors respectively, and Mj is the fermion mass. The functions a^, a R , b L , 
b R , cl and c R are the chiral projections of the Lorentz-invariant functions a, b, c 
defined in the following way: 

a = a^L + a R R, (2. 2) 

b = b L L + b R R, (2. 3) 
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C = ClL + CrR. 

These functions depend on the Lorentz scalars u and k defined by u 
and k = [{K ■ u) 2 - K 2 ] 1 / 2 . 

The inverse fermion propagator is given by 



S-\K) = (L + fR- A L L - A R R, 



where: 



C? = (1 + a L )K» + b L u\ 
W = (1 + a R )K» + b R u^, 
A L = M r (l-c L ), 
A R = M I (l-c R ). 
The fermion propagator follows from the inversion of (2. 5) 

1 



(2. 4) 
(K-u) 



(2. 5) 

(2. 6) 
(2. 7) 
(2. 8) 
(2. 9) 



(£ 2 f - A L A R () L + (K 2 ^ - A L A R $) R + ArL(£-U- A l A r + - [£ 

+ A l r(jC-&- A l A r + X - [fi, R 



The poles of the propagator correspond to values of u and k for which the 
determinant D in (2. 10) vanishes Pf, 



D(u, k) = C 2 W - 2A L A R (C ■ K) + (AlAr) 



0. 



(2. 11) 



L 

(2. 10) 



In the rest frame of the plasma u = (1,0), Eq. (2. 11) leads to the fermion dis- 
persion relation for a gauge theory with chiral and parity violations (for instance, 
the MSM after the ESSB): 

{[w(l + a L ) + b L f -k 2 [l + a L ] 2 } + a R ) + b R f - k 2 [1 + a R } 2 } - 
2M 2 (1 - c L )(l - c fl ) {(^ 2 - P)(l + a L ) (1 + 0ji ) + uj[(1 + a L )b R + (1 + a R )b L ] + 6 L ^} 



+M / 4 (l-c L ) 2 (l-c R ) i 



0. 



(2. 12) 



Equation (2. 12) describes in a coupled way the behaviour of the left and right 
components. This equation is of fourth order in u, so that for each value of 
k there are four solutions for u. Therefore four branches (two normal and two 
abnormal) in the phase space are obtained from the dispersion relation. It is 
not possible to associate a specific chirality with the quasi-particles because the 
chiral symmetry is broken. 

For massless fermions, the chiral symmetry is preserved by the theory. The 
functions and A# in (2. 11) are equal to zero and the poles of the propagator 
are found from £ 2 3? 2 = 0. Thus, for a chirally invariant theory with parity 
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violation, as is the case of the MSM before the ESSB, the dispersion relations are 
given by 

+ a L ) + b L } 2 - k 2 [l + a L } 2 = 0, (2. 13) 

[u(l + a R ) + b R ] 2 -k 2 [l + a R f = 0. (2. 14) 

Left- and right-handed components now obey decoupled relations. Equations 
(2. 13) and (2. 14) are of second order in uj; as a consequence, for each mo- 
mentum value k there are two values of u re-solving each equation. For each 
fermion chirality, two branches (one normal and one abnormal) are solutions of 
the dispersion relation. Four branches describe the fermionic quasi-particles of 
the plasma. The quasi-particles have a specific chirality, in contrast with what 
happens in the chiral violation case. 

For a parity-symmetric gauge theory, the projections of the Lorentz-invariant 
functions over the left- and right-handed components coincide. Setting ai = 
a R = a, bi = b R = b and cl = c R = c in (2. 12), the dispersion relation for 
a parity-invariant gauge theory with chirality violation (for instance, QCD with 
massive quarks) is obtained: 

HI + a) + b} 2 -k 2 [l + a] 2 - M 2 [l - c] 2 = 0. (2. 15) 

The left- and right-handed components of the dispersion relation are super- 
imposed in this case, and the fermionic excitations are only described by two 
branches. 

If chiral symmetry is preserved by a parity-invariant gauge theory, as is the 
case of QCD with massless quarks, the dispersion relation is given by (2. 15), 
with M 7 = 0, 

[w(l + a) + b] 2 - k 2 [l + a} 2 = 0. (2. 16) 

3 Real part of the fermion self-energy 

The finite temperature contribution of the one-loop fermion self-energy is cal- 
culated in this section. The MSM Feynman diagrams that contribute to the 
one-loop fermion self-energy have the form of the generic diagrams of Figs. 1 
and 2. These generic diagrams are calculated, in the Feynman gauge, using the 
real time formalism of thermal field theory |14| . It is straightforward to see that 
tadpole diagrams represented by Fig. 1 only contribute to the imaginary part 
of the self-energy, as is shown in Appendix A. The diagrams that contribute to 
the real part of the self-energy are represented by the generic diagrams of Fig. 2. 
To calculate the poles of the fermion propagator, we are only interested in the 
contribution to the real part of the fermion self-energy from these diagrams. The 
same calculation is also performed in the Landau gauge (Appendix B). 
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Finite temperature Feynman rules for vertices are the same as those in the 
MSM at T = 0, while the thermal propagators for massive fermions S(p), massive 
gauge bosons D^ u (p) and massive scalar bosons D{p) are given by 



S(p) = + M x 



1 



9nu 



D(p) 



p 2 - Mf + ie 
1 

p 2 - M$ + ie 
1 



+ iT f (p) 

- «r b (p) 



- iT 6 (p), 



(3. i; 



(3. 2) 



(3. 3) 



p 2 - Ml + ie 

where the T dependence is introduced through the functions Tf{p) and r&(p) 
defined by 



T f (p) = 2tt5(p 2 - Mf)n F (p 



u 



(3. 4) 
(3. 5) 



T b (p) = 2n6(p 2 -M 2 )n B (p-u), 
where np(p ■ u) and n B (p ■ u) have the following form: 

n F (p ■ u) = [e [ l p - u| / T] + I]" 1 , (3.6) 

n B (p-u) = [e [ ^ T] - l]" 1 . (3. 7) 

The contribution to the self-energy from the generic gauge boson diagram of 



Fig. 2a is 



(27T) 



(3. 



where g is the coupling constant and C c is the quadratic Casimir invariant of 
the re- presentation. Inserting expressions (3. 1) and (3. 2) into (3. 8), it is 
possible to write: S(i^) = E°(i^) + Y!(K), where YP(K) is the zero temperature 
contribution and S'(iT) is the finite temperature contribution. Keeping only the 
real part ReE'(iT) of the finite temperature contribution, we obtain 



ReZ'(K) = 2g 2 C c J ^-{j>^-2M_ 



(27T) 



T 6 (P) 



(p + iT) 2 - Mj 2 



Ml 



, (3. 9) 



where the denominators are defined by their principal value. Applying Weldon's 
procedure and after performing the integrations over p and the two angular 
variables, the following novel integrals over the modulus of the three-momentum 
p = \p\ are obtained: 



-Tr(ReS') 
4 ' 



-2g 2 C c M\ 



dp p ' 1 1 

— 27 — L t{p)nB{e 2 ) + —Lt(p)n F (ei) 
2 k Le 2 ei 



in 



(3. 10) 
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1 r c 
-Tr(^ReE') = g 2 C c 
4 Jo 



+ — 

ei 

dp p 



1_ f (u 2 -k 2 - A) 
i V 2 
1 /(u; 2 -£; 2 -A) 



L+(p) +Apk \n B (e 2 ) 



L+(p)+4pk)n F (e 1 ) 



(3. 11) 



8vr 2 jfe 

with the logarithmic functions given by 



L 2 (p) + — L t (p) n B (e 2 ) + L 1 (p)n F (ei) 

(3. 12) 



Lf (p) = log 



L 2 (p) = log 



cu 2 - k 2 - A - 2 ei cj - 2A;p' 
u 2 - k 2 - A - 2eiu; + 2£;p 

cj 2 - /c 2 + A + 2e 2 cj + 2£;p' 
cu 2 - A; 2 + A + 2e 2 cj - 2£;p 



±log 



± log 



uj 2 - k 2 - A + 2eicj - 2fcp~ 

cu 2 - k 2 - A + 2eiu + 2kp 

(3. 13) 
lu 2 - k 2 + A - 2e 2 u; + 2£;p 
cj 2 - k 2 + A - 2e 2 c<; - 2kp 

(3. 14) 

where k = \k\, e x = (p 2 + M 2 ) 1 / 2 , e 2 = ( p 2 + M 2 ) 1 / 2 , an d A = Mf - M 2 . These 
expressions for Mi = M 2 = are identical to those derived by Weldon |J and 
for Mi ^ and M 2 = to those derived by Petitgirard ||. 

The contribution to the fermion self-energy from the generic scalar boson di- 
agram shown in Fig. 2b has the same form as the gauge boson contribution. The 
factor 2g 2 C c in (3. 9) is replaced by \l 2 \C' LR , where I is the Yukawa coupling 
constant and C' L R are given in terms of the matrices of Clebsch-Gordan coeffi- 
cients P. The generic scalar boson contribution for left and right chiralities is 



proportional to (PP 



l" in in' an d (P P 



C'r.S. 



C' R 5 nn i respectively 



The integrals (3. 10)- (3. 12) are related to (2. 1) through 



^Tr(fiRe £') - w~Tr(^ Re Tf) 



2 - fc 2 )~Tr(^ReE') - u~Tr(fiReTf) 
c(u,k) = --^-^Tr(ReE'), 



(3. 15) 
(3. 16) 
(3. 17) 



where the Lorentz-invariant functions a (a;, k), b(ui, k) and c(u, k) can be written 
in terms of their chiral projections, as was shown in (2. 2)-(2. 4), and Mj is the 
mass of the external fermion. 

It is convenient to rewrite (3. 15)-(3. 17) in the following form: 



a{uj,k) = g 2 C c A(M u M 2 ) 
b(u 1 k)=g 2 C c B(M 1 ,M 2 ) 



(3. 18) 
(3. 19) 
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Mi 

c{u,k)=2g 2 C cl ±C{M 1 ,M 2 ) 



(3. 20) 



with 

A(M U M 2 ) 



1 f°° dp 



k< 



I 

Jo 



8tt 2 
+ 



2fc 7 2 2 {p) - T 2 (p) + 77 



{u> 2 -k 2 - A) p up V 

2l ^ (p) "T Ll(p) + - 



rip^)) (3. 21) 



B(M 1 ,M 2/ 



1 Z" 00 dp 



A; 2 Jo 8tt 2 



feu 2 - A; 5 



A: 



-pL 2 (p) + 



(uj 2 -k 2 + A)u p 



k e 2 



+ 



(oj 2 - k 2 
k 



-pL l (p) - 



(uo 2 - k 2 - A) uj p 
2 ~k~ x 



Ltip) ~ 



Aujp 2 



Atop" 2 
n F (e 1 ) 



) , (3- 22) 



C(M 1; M 2 ) = 11°°^ (^Lt(p)n B (e 2 ) + ^L+(p)n F (ei) ) . 



., I <-• ■> > ( ' i ' ■' I > ■ ! l-i I \ /' / " /• \ > I . I I ■ (d. 26) 

on 2 \e 2 ei J 

Since each one-loop diagram contribution to the real part of the fermion self- 
energy is proportional to the generic gauge boson contribution (3. 9), it is possible 
to write 

a(u,k) = J2fnA(M u M 2 ), (3. 24) 



b(u,k) = J2f n B(M 1 ,M 2 ), 

n 

c(u,k) = J2g n C(M l ,M 2 ), 



(3. 25) 
(3. 26) 



where the Lorentz-invariant functions are obtained by adding up the contributions 
n coming from each diagram involved in the calculation; M 1 and M 2 refer to the 
internal fermion and to the boson mass of the loop diagram, respectively. 



4 Fermionic dispersion relations 

In this section, the one-loop dispersion relations at finite temperature for 
quarks, charged leptons and neutrinos in the MSM are derived. The chiral pro- 
jections of the Lorentz-invariant functions in the broken and in the unbroken 
phases are calculated following the procedure explained in the last section. The 
unbroken phase Lorentz-invariant functions calculated in the Landau gauge are 
presented in Appendix C. 

In the MSM, the W and Z gauge bosons and the fermions acquire mass after 
the ESSB. These masses are written as functions of v in the form: = g 2 v 2 /A, 
M z — {g 2 + 9 >2 )v 2 / A and Mp = lpv 2 /2, where g and g' are the gauge coupling con- 
stants of SU(2) L and U(1) Y respectively and l F is the Yukawa coupling constant 



9 



of the fermion F to the Higgs boson. The mass of the Higgs boson is M\ = pv 2 /2, 
where p is a free parameter of the Higgs potential. The thermal dependence of 
the boson and fermion masses is given through v by v(T) ~ v (l - T 2 /T c ) 1 ' 2 , 
where Tc ~ 100 GeV is the critical temperature of the electroweak phase tran- 
sition and vq = 246 GeV is the vacuum expectation value of the Higgs field at 
T = 0. The latter expression is valid for values of T near and below Tc- 



4.1 Quark dispersion relations 

As we have already mentioned, the damping rate of the quasi-particles in the 
plasma is proportional to the imaginary part of the fermion self-energy. The 
damping rate is induced by the incoherent scattering of the quasi-particles with 
the thermal medium. The scattering of the quasi-particles associated to quarks 
is generated by the strong and electroweak thermal interactions. The most im- 
portant contribution to the damping rate comes from QCD thermal interactions. 
It has been calculated at leading order for zero quasi-particle momentum in JIT 



giving: •Jqcd ~ 0.15g 2 T, i.e. ~ 19 GeV at T = 100 GeV, where g s is the strong 
coupling constant. 

Assuming that, in the vicinity of k = 0, the damping rate is not changed 
appreciably by the electroweak thermal interactions, the finite lifetime of the 
quasi-particles is ~ 1/40 GeV -1 . Since the group velocity is ~ 1/3, the mean free 
path is ~ 1/120 GeV -1 . Thus, the quantum spatial coherence is lost for points 
separated by a distance > 1/120 GeV -1 . We show below that the electroweak 
corrections to the branches for light quarks are of the order of the quark masses 
Mj. Thus, the imaginary part of the QCD self-energy is much larger than the 
electroweak contributions to the real part of the self-energy for all quarks but 
the top. This means that coherent electroweak quantum phenomena associated 
to quarks (for instance, CP violation) cannot take place in the plasma because 
the lifetime of the quasi-particles is very small with respect to the characteristic 



time of the electroweak processes 0, [12|]. When a quasi-particle scatters, the 
electroweak coherence is lost. 



4.1.1 Broken phase 

The thermal background interacting with quarks in the plasma is constituted 
by massive quarks of three flavours with colour charge, massive Higgs bosons, 
massive W and Z gauge bosons, photons and gluons. Quarks of electric charge 
2/3 are called up quarks and those of electric charge —1/3 down quarks. The 
one-loop diagrams with photon, gluon, Z gauge boson, neutral Goldstone and 
Higgs boson, which contribute to the quark self-energy, do not induce a change 
of flavour in the quark. The diagrams with an exchange of charged Goldstone and 
W gauge bosons induce a flavour change in the incoming quark / to a different 
outgoing quark F. In the latter contributions, the flavour i of the internal quark 
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(inside the loop) runs over the up or down quark flavours according to the type 
of the external quark (outside the loop). The one-loop contribution from gauge 
boson and scalar boson to the real part of the self-energy gives 

a L (u,k) IF = [fA(M I ,0)+f!A(M I ,M z )+f H A(M I ,M H )]5 IF +Y,f^A(M lJ M w ) J 

R R • R 

(4. 1) 

b L (uj,k) IF = [fB(M I ,0)+f!B(M I ,M z )+f H B(M I} M H )]5 IF +Y,f^B(M l} M w ), 

R R i R 

(4. 2) 

c L (u,k) IF = [2fC(M I ,0)+g z C(M Ij M z )-f H C(M I ,M H )]5 IF +Y,9fC(M l ,M w ). 

R I R 

(4. 3) 

The integrals (3. 21)-(3. 23) will be evaluated with the appropriate values of M 1 
and M.2 for each contribution. The coefficients f n and g n are: 

/ = \gl + g 2 Q 2 isi, (4- 4) 

ft = C(T? ~ Qtsl? + fi = 9 ^Q]si + (4. 5) 

C w o C w o 

f H = (4- 6) 

fY = y(l + %)Kt<KiF, f% = 9 -^Kf t K lF , (4. 7) 

9 Z = - 2 -^(Tf-Q I sl)Q I sl + 9 -^, (4.8) 
c w o 

9Y = 9 -^Kf t K lF , gj = 9 ^Kt lKlF) (4. 9) 

where K represents the Cabibbo-Kobayashi-Maskawa (CKM) matrix, A/ is 
defined as A/ = Mj/Mw, g is the weak coupling constant, Qi is the quark charge 
and s w (c w ) is the sine (cosine) of the electroweak mixing angle 6 W . The Yukawa 
coupling constant lj is related to A/ through lj = g\j /y/2. The substitution of 
the expressions (4. l)-(4. 3) in (2. 12) results in a non-diagonal matrix because of 
the presence of mixing in the CKM matrix. From the quark dispersion relation 
obtained, the branches are computed. The mass non- degeneracy of the quark 
spectrum and the different Yukawa coupling constants for the distinct quark 
flavours introduce the flavour non- degeneracy of the branches. The branches 
associated to each quark flavour are different because of the presence of Mi, A/, 
Qi and the CKM matrix in (4. l)-(4. 3). 

The spectrum of quasi-particles for massive quarks with only QCD interac- 
tions was studied in || [8|. The quasi-particles are described by one normal and 
one abnormal branch. If electroweak effects are introduced, as they are here, the 
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mentioned spectrum is modified in a complicated way. Although the electroweak 
corrections are quantitatively small for all flavours but the top, some important 
consequences are derived from these effects. In particularly, as was explained 
in section 2, four branches are present (two normal and two abnormal) instead 
of the two QCD ones. The four branches for the s and c quarks are shown in 
Figs. 3 and 4. Because the chiral symmetry of the theory with massive quarks is 
broken, it is impossible to associate a specific chirality to the branches. The top 
quark acquires a large mass during the electroweak phase transition. We note a 
particular behaviour of the top quark branches with the T variation: a couple of 
branches lowers down appreciably as T is lowered (M t increases) as is shown in 
Fig. 5, by means of the T variation of the thermal effective masses. To try to 
understand this special behaviour of the top quark branches, we will first study 
the much simpler QCD case. 

It is very easy to check that (4. l)-(4. 3) reduce to QCD functions when g is 
fixed to zero in the coefficients (4. 4)-(4. 9). That is for = a R = a, bi J = bn = b, 
cl = c r = c, 

a(u,k) = f'A(M I ,0), (4. 10) 



b(u, k) 
c(lu, k) - 



= /'B(Afj,0), 
2/'C(M /s 0), 



(4. 11) 
(4. 12) 



where /' = |g 2 . Inserting these expressions in (2. 15), the dispersion relation for 
massive quarks in QCD is obtained. We have checked the validity of the analytic 
approximations performed in || for light massive quarks (Mj << T). In this 
case, the analytic forms of the normal and abnormal branches for k «T are: 
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M+ + 



1 



M 



Lead 



Mlead 

3M, 



+ 



2Mr 



3M 2 



M 2 Lead 



3M 4 



k 2 + 0(k 4 ), (4. 13) 



M_ + 



1 



M 2 Lead + Mi 



Mlead 

3M_ 



'3M 2 



MLadX 



2Mj 



3M_ 



k 2 + 0(k 4 ), (4. 14) 



respectively, with M+ = [(Mf +4M| ead ) 1 / 2 +M / ]/2 and M_ = [{M 2 +AM 2 ead ) 1 ^ 2 - 
M/]/2, where M 2 Lead = f'T 2 /8. If k = then uo = M + for the normal branch 
and uj = M_ for the abnormal branch. These different values of the energy are 
interpreted as thermal effective masses. For values of k close to zero, the normal 
branch is above (M 2 + 4M| ead ) 1 / 2 /2 by Mj/2 while the abnormal one is below 
that value by the same amount. The separation between the branches is Mj 
and this increases with the quark mass. When the electroweak corrections are 
taken into account, a similar effect is also present. An energy gap of order ~ Mj 
appears between the minimum and the maximum of the hyperbolic branches 
(see Figs. 3 and 4). If the quark mass is large, the energy gap is larger than 
the separation between the thermal effective masses. As a consequence, the 
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normal and abnormal branches with the same thermal effective mass appear 
closer. For some T, the thermal effective mass of smaller energy can be negative, 
as is suggested by Fig. 5. 

More works about the dispersion relation for massive quarks in QCD can 
be found in the literature. Analytic expressions for the dispersion relation were 
calculated, in the limit Mf + k 2 » Ml ead , by Seibert 0. In the limit Mf >> 
M\ ead as was done in |16| , the approximation for the dispersion relation oj 2 = 
k 2 + Mf + Mf ead is valid. The properties of the quark spectrum were studied by 
Lebedev and Smilga [17j considering a dynamical thermal mass for the quarks in 
the tree-level matrix Green's Functions. 



4.1.2 Unbroken phase 

With the masses of the particles equal to zero, the expressions (4. l)-(4. 3) 
are reduced to the chiral projections of the Lorentz-invariant functions in the 
unbroken phase: 

a L (u,k) IF = l IF A(0,0), (4. 15) 

b L (u,k) JF = l IF B(0,Q), (4. 16) 

a R {u,k) IF = r IF A(0,0), (4. 17) 

b R {u,k) IF = r IF B{0,0), (4.18) 
where the coefficients h F and rj F are given by 



IF 



1 + 



tan 2 6» 



IV 



27 



$IF + 



K 



it 



X-K iF 



riF = S 



I 3 + ~Y~ 



Q^tan 2 6 W + -Xj 



6 



IF- 



(4. 19) 
(4. 20) 



The branches for massless quarks are computed from the dispersion relation ob- 
tained by substituting (4. 15)-(4. 18) in (2. 13) and (2. 14). Two normal branches 
(one left-handed and one right-handed) and two abnormal branches (one left- 
handed and one right-handed) describe quasi-quarks and quasi-holes respectively. 
The different Yukawa coupling constants for the distinct quark flavours introduce 
the flavour non-degeneracy of the branches. In Figs. 3 and 4, the s and c quark 
branches are presented. 

Calculating for u,k T the leading 0{T 2 ) contribution to the expressions 
(4. 15)- (4. 18), it is possible to observe analytically the chiral non- degeneracy of 
the quasi-particle spectrum. The integrals A(0, 0) and -8(0, 0) in this limit allow 
us to write 

/ 7 \ i 1>ifT r 

a L (u, k) IF = - — 1-F 



8 k 2 



b L (u,k) 



HifT 



IF 



8 k 



OJ 

2k 



y w k' \2kJ 



(4. 21) 



(4. 22) 



13 



a R (u, k) 



lr IF T 



2 r 



IF 



8 k 2 



UN 



2kJ 



b R (u, k)i F 
where F(x) is defined by: 



lr IF T 2 



CO 

T 



Fix) 



x 
2 



log 



H (- 
w 

x + 1 
x — 1 



to, 



(4. 23) 
(4. 24) 

(4. 25) 



In this limit, the quark self-energy has been calculated for QCD in [TS j and 
generalized for the MSM in HJ. The following analytical expressions are derived 
by replacing (4. 21)-(4. 24) in (2. 13) and (2. 14): 



wlr - k 



Ml R 



1 + 



1 



k 



Ilo I UL > R + k ' 
2 V^-B ~ 



(4. 26) 



where we have defined M? 



, L - l u T 2 /8 and M 2 R = r u T 2 /8. The normal and 
abnormal branches for small momentum [k « Ml,r) are given by 



k 2 



k k 2 

UJlA k ) = M L,R ~ q + 



+ C(A; 3 



(4. 27) 



(4. 28) 



respectively. For k = 0, cl^b = and then Ml,r are interpreted as thermal 

effective masses. For each chirality, one normal and one abnormal branch describe 
the quasi-particle spectrum. For high momentum (k >> M^), the branches are 
approximated by 



^L,R{k) = k + 



Ml R 



Kr 
2fc 3 



log 



' 2k 2 ' 
MIr, 



+ 



UL,R{k) = k + 2k ■ exp 



2k 2 



+ 



(4. 29) 



(4. 30) 



For k » Ml^r the branches become the ordinary fermionic branch. The leading 
0(g 2 T 2 ) expressions (4. 21)-(4. 24) are gauge independent and they give place 
to the gauge invariant dispersion relations (4. 27)-(4. 30) ||. Supposing that 
the propagation of quasi-particles is along the z axis, it is possible to see that 
their group velocity is of ±1/3. For k z and a given flavour, quasi-particles of 
a specific chirality have two possible helicities. The group velocity is given by 
the eigenvalue of ~ l/37 5 <7 2 = l/3'j 5 a z (k z ) 2 , i.e. l/3xhk z , where h is twice the 
helicity, x the chirality and k z = k z /\k z \ 0. Consequently, the group velocity 
for quasi-fermions has the same sign as the momentum when chirality equals 
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helicity. For quasi-holes, the group velocity and the momentum have opposite 
signs. These fermionic states have flipped chirality/helicity and at zero T they 
have negative energy. They become physical at finite T due to a shift in energy 
of the branch. 

Turning off the electroweak interactions in the Lorentz-invariant functions 
(4. 15)-(4. 18), the QCD case for massless quarks is recovered, with 

a(u,k) = f'A(0,0), (4. 31) 

b(u,k) = f 5(0,0). (4. 32) 

The corresponding dispersion relations follow from the substitution of (4. 31) and 
(4. 32) in (2. 16). In the limit u, k << T the branches are given by (4. 27) and 
(4. 28), with u L = u R = u and M L = M R = M Lead in accordance with p[ p|, pp. 



4.2 Charged- lepton dispersion relations 

The scattering between the quasi-particles associated to charged leptons and 
the thermal medium is governed by the weak and electromagnetic thermal inter- 
actions. It is expected for k = 0, in analogy with QCD calculations ||11|| , that 



the damping rate at leading order originated only from weak or electromagnetic 
thermal interactions were 'Jew ~ 9 2 T. The weak and electromagnetic leading 
contributions of the real part of the self-energy are both of order gT. This means 
that the energy of the quasi-particles, of order gT, is large in comparison with 
the damping rate, i.e. it contains one additional power of g. Thus, it is possible 
to speak of the existence of quasi-particles as coherent excited states. Coherent 
electroweak quantum phenomena associated to charged leptons can take place 
in the plasma because the lifetime of the quasi-particles is of the order of the 
characteristic time of the electroweak processes. 



4.2.1 Broken phase 

The thermal background interacting with charged leptons in the plasma is 
constituted by massive charged leptons of flavour e, massless neutrinos of flavour 
e, massive Higgs bosons, massive W and Z gauge bosons and massless photons. 
The diagrams with charged Goldstone and W bosons change from the external 
charged lepton to the associated neutrino flavour inside the loop. The one-loop 
contribution from scalar and gauge boson sectors to the real part of self-energy 
allow us to write 

a L (u,k) = dA(M e ,0) + d z L A(M e ,M z ) + d H A(M e ,M H ) + df A(0,M w ), (4. 33) 

R R R 

b L (u,k) = dB(M e ,0) + d z L B(M e ,M z ) + d H B(M e ,M H ) + d^B(0,M w ), (4. 34) 

R R R 

c l (lu, k) = c r (lu, k) = 2dC(M e , 0) + y z C(M e , M z ) - d H C(M e , M H ), (4. 35) 
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with the coefficients d and y given by 



d = g 2 s 2 w} (4. 36) 



n 2 1 n 2 \ 2 n 2 n 2 \ 2 

4 = 9 ~r(\ - si) 2 + g -j±, dl= 9 -si + (i. 37) 



w * w 

,2 \2 



d H = £A (4. 38) 



# = £ < = ^ (4- 39) 



f = -^-4K + V £ . (4- 40) 



where A e = M e /M w . Substituting the functions (4. 33)-(4. 35) in (2. 12), the 
branches are computed from the obtained dispersion relation. Because the chiral 
symmetry in this case is broken, it is impossible to associate a specific chirality 
to the branches. Flavour non-degeneracy exists since the mass non-degeneracy is 
present in (4. 33)- (4. 35). 

The spectrum of quasi-particles for charged leptons in QED was studied in 



HJ. For this case, the quasi-particles are described by one normal and one abnor- 
mal branch, and left- and right-handed branches are superposed. The breaking 
of the parity symmetry by the weak interactions causes the doubling of the QED 
branches. Thus, for charged leptons in the MSM, the quasi-particles are described 
by two normal and two abnormal branches, as is shown in Fig. 6 for the r lepton. 
As in the quark case in the broken phase, an energy gap between the hyperbolic 
branches appears. 

Following a procedure analogous to the one performed for massive quarks in 
QCD, the Lorentz-invariant functions for massive electrons in the QED case are 
obtained by turning off the weak interactions into the functions (4. 33)- (4. 35). 
The expressions derived are similar to expressions (4. 10)- (4. 12), with /' changed 
to d and Mj to M e . The analytic form of the two branches, for k « T, is given 
by (4. 13) and (4. 14). 

Some aspects of the dispersion relation for massive electrons in QED have 
been studied at finite temperature, considering an effective chemical potential in 



21 1 . The spectrum of quasi-particles in a hot relativistic electron plasma was 



studied in p2|, where it is shown numerically that the full spectral weight is not 



gauge independent. By comparing one- loop computations in the Feynman and 



in the Coulomb gauges, it was shown in |22| that the location of well-defined 
poles in the propagator is gauge independent at leading order since the poles 
correspond to physical excitations, while the widths of the quasi-particles are 
gauge dependent. The dispersion relation at zero temperature has been studied 
23|. 



in 
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4.2.2 Unbroken phase 

Expressions (4. 33)-(4. 35) with particle masses equal to zero reduce to the 
chiral projections of the Lorentz-invariant functions in the unbroken phase: 

a L (u,k) = tA(0,0), (4. 41) 

b L (u,k) = tB(0,0), (4. 42) 

a R (u,k) = uA(0,0), (4. 43) 

b R (u,k)=uB(0,0), (4. 44) 



where the coefficients t and u are 

1 



t = TCCtw 



u = Ana 



w 



+ -r + A 

c 



A 2 

tan 2 Ow + y 



(4. 45) 



(4. 46) 



The branches for massless charged leptons are computed from the dispersion re- 
lation derived by substituting (4. 41)-(4. 44) in (2. 13) and (2. 14). Two normal 
branches (one for each chirality) and two abnormal branches (one for each chiral- 
ity) describe quasi-particles and quasi-holes respectively (see Fig. 6). The quasi- 
particle spectrum is flavour non-degenerate, but, because A 2 in (4. 45) and (4. 46) 
is so small for all charged lepton flavours, the differences between flavours are very 
small. The type of quasi-particle characterized by the abnormal branch has been 
called the "plasmino" by Braaten and Pisarski P^fl . In the limit u, k <C T, the 
analytical forms obtained for the dispersion relation are similar to those obtained 
for massless quarks. For k « M^^, the branches are given by (4. 27) and (4. 28) 
with M 2 L = tT 2 /8 and M 2 = uT 2 /8. 

4.3 Neutrino dispersion relations 

The damping rate of quasi-particles associated to neutrinos is induced by the 
weak thermal interactions. In analogy with the assumption in the charged- lepton 
case, the damping rate (leading order) of neutrino quasi-particles for k = is 
7evk ~ g 2 T. The weak leading contribution of the real part of the self-energy 
is of order gT. The energy of the quasi-particles, of order gT, compared to 
the damping rate is large. The quasi-particles exist as coherent excited states. 
Coherent weak quantum phenomena associated with neutrinos can take place 
in the plasma because the lifetime of the quasi-particles is of the order of the 
characteristic time of the weak processes. 
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4.3.1 Broken phase 



The thermal background interacting with neutrinos in the plasma is consti- 
tuted by massless neutrinos of flavour e, massive charged leptons of flavour e 
and massive W and Z gauge bosons. The one-loop contribution from scalar and 
gauge boson sectors to the real part of the self-energy lead to 

a L (uj, k) = v z A(0, M z ) + v w A(M e , M w ), (4. 47) 

b L (cu, k) = v z B{0, M z ) + v w B{M e , M w ), (4. 48) 

with the coefficients v z and v w given by 

vi = (4. 49) 

c 

v w = na w [2 + X 2 e ], (4. 50) 

where A e = M e /Mw- The projection of the Lorentz-invariant functions is only 
over the left-handed component. Substituting (4. 47) and (4. 48) in (2. 13), the 
branches are computed from the dispersion relation obtained. A normal and an 
abnormal branches of left-handed chirality describe the quasi-particles as is shown 
in Fig. 7. Flavour non- degeneracy exists from the mass non-degeneracy present 
in (4. 47) and (4. 48), for different neutrino flavours. Because M e << Mw in 
the integrals of (4. 47) and (4. 48), the differences for different flavours are very 
small. 

The neutrino dispersion relation presented here was studied by D'Olivo, Nieves 
and To-rres in [^J for temperatures of the plasma much lower than the W and 
Z masses. This approximation allows to neglect the thermal dependence for the 
gauge boson propagators. These authors performed the calculation in a general 
gauge and they found that the neutrino self-energy depends on the gauge pa- 
rameter, but the dispersion relation is independent at leading order in g 2 /M w . 
Several works, previous to G5j, on this subject can be found in the literature. In 



27], the same calculation was performed but using the unitary gauge. The 



computation of the neutrino thermal self-energy was considered in PR . 



4.3.2 Unbroken phase 

Expressions (4. 47) and (4. 48) with gauge boson masses equal to zero reduce 
to the left-handed projection of the Lorentz-invariant functions in the unbroken 
phase: 

a L (u,k) = L4(0,0), (4. 51) 

b L (u,k) = 1B(0,0), (4. 52) 

with the coefficient t given by (4. 45). Substituting (4. 51) and (4. 52) in (2. 13), 
the branches are computed from the dispersion relation obtained. A normal 
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and an abnormal branch of left-handed chirality describe the quasi-particles (see 
Fig. 7). The quasi-particle spectrum is flavour non-degenerate, but, since \ 2 e is 
very small for all neutrino flavours, the differences for the different flavours are 
very small. The variation with the temperature of the thermal effective masses 
for both phases is presented in Fig. 8. 

The branches for left-handed neutrinos are the same as those for left-handed 
charged leptons. This result reflects the presence of the SU(2) l symmetry in the 
MSM before the ESSB. The neutrino dispersion relation in the unbroken phase 



was considered with a non-zero chemical potential in |30(]. The same authors in 



a later work |]31| studied the dispersion relation in the case where the charged 
leptons have a running mass obtained using the renormalization group equations. 



5 Remarks about gauge invar iance 

The gauge invariance of the one-loop dispersion relations presented in this 
work is studied in this section. The numerical results for the dispersion relations 
have been obtained considering all the terms in the Lorentz-invariant functions 
(3. 18)-(3. 20). This fact implies the introduction of a small gauge dependence 
originated from the non- leading order terms of the fermion self-energy |H| . These, 
in turn, lead to sub-leading contributions to the dispersion relations. Indeed, by 
confronting QCD computations in the Feynman and Landau gauges, we show that 
the sub-leading contributions to the dispersion relations are gauge dependent. We 
also study the effect of these contributions in our numerical results. 

In Appendix B, we calculate the finite temperature contribution to the one- 
loop real part of the fermion self-energy in the Landau gauge ReS'^ (K). We 
show in (B. 4) that it is possible to write: ReE' £ (K) = ReE'(K) + Re5T,'(K), 
where ReH'(K) is the real part of the thermal contribution to the self-energy in 
the Feynman gauge (3. 8) and ReSH'(K) is the real part of the gauge dependent 
thermal contribution of ReH' £ (K) (B. 5). The Ward identities ensure that the 
gauge dependent contribution disappears if the external leg is projected on the 
zeroth order pole ||32|| . This can easily be seen because Re8T,'(K) can be written 
as proportional to I/C— M q . In QCD at finite temperature, for instance, the shift 
of the quark propagator pole position with respect to the pole position at zero 
temperature is of order g s T. This means that the argument of the integral in 
(B. 5) is proportional to g s T, for this case, because this integral is proportional 
to $T— M q . Consequently, the Re5Y,'(K) contribution to the dispersion relation 
is of a higher order in g s with respect to the ReH'(K) contribution. 

For example, in the QCD case for massless quarks, Re8Ti'(K) vanishes on 
the zeroth order pole corresponding to uj = 0, for k = 0. So, ReS'E'(K) con- 
tribution is a non- leading order in uo/T. These non- leading order contributions 
lead to a sub-leading contributions to the thermal effective mass. This fact can 
be guessed because the leg projection on the first order pole means uj/T ~ g s . 
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This conclusion is in accordance with |32" |, where it is formally shown that non- 
leading contributions to the quark self-energy allow gauge dependent sub-leading 
modifications to the dispersion relation. 

Let us now show in an explicit calculation that the non-leading contributions 
to the one-loop QCD self-energy are sub-leading in the dispersion relations. By 
comparing one-loop computations in the Feynman and in the Landau gauge for 
massless and for massive quarks, we will show that indeed the sub-leading con- 
tributions to the dispersion relations are gauge dependent. First, we consider 
massless quarks in QCD. The Lorentz-invariant functions in the Feynman gauge 
(3. 18) and (3. 19) for this case are 



a Q (u,k) 



M 2 Lead 



UJ 

2k- 



+ 



4! 



b Q (uj,k) 



M 2 Lead 
k 



uj k uj / uj \ 
k + [ -u~V \2k) 



u 2 h + k 2 h - 2ujI 3 

2 -k 2 



+ 2A 



,UJ^ 



k a 



a 2 



(5. 1) 



(5. 2) 



where M\ ead = g 2 T 2 /6, A 2 = g 2 /12n 2 and g s stands for the strong coupling 
constant. The function F(x) is defined by (4. 25) and the integrals I 1; I 2 , -^3 are 
given by 

POO 

h(uJ,k,T) = / dp Lf{p) [n B (p)+n F (p)} , (5. 3) 



I 2 (uJ,k,T) 
I 3 (uj,k,T) -- 



dp Lj(p) [n B (p) -n F (p)\, 
dp pLf(p) [n B {p) + n F { 



(5. 4) 
(5. 5) 



where ns(p) and n F (p) are defined by (3. 6) and (3. 7) respectively. The loga- 
rithmic functions Lf(p) are 



Lf{p) = log 



2p - (u - k) 
2p-(u + k) 



±log 



2p + (uj + k) 
2p + (u) - k) 



(5. 6) 



Substituting the expressions (5. 1) and (5. 2) into (2. 16), it follows that 



Kead 



k 



uj — k 



- log 

2 6 



uj + k 
uj — k 



A 2 (uh - kl 2 - 2/ 3 ) . 



(5. 7) 



For k — > 0, the integrals Ji, J 2 , J 3 can be written as Ii(uj,k,T) = kFi(u,T) + 
0(k 3 ), where i = 1,2, 3. The reason of this is that (5. 3)-(5. 5) are antisymmetric 
for the change k — > —k. It was not possible for us to obtain analytic expressions 
for the integrals (5. 3)-(5. 5). They were computed numerically for different values 
of uj, k and T, in the limit k — > 0. By means of a fit of the numerical results 
obtained, we have found an analytical form for the integrals. Numerically we 
have proved that (uil\ — 2/3) = 0(k 3 ), i.e. (ujFi — 2F 3 ) = 0. For this reason 
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(ojI\ — 2/3) vanishes in (5. 7), when we take the limit k << uj, T > uj, divide by 
k 2 and set k = 0. Then, we obtain the following relation: 



uj 



A 2 F 2 (uj,T) 



M 2 

lvl Leadi 



(5- 



wheref] 



F 2 (u;,T) 



iO 



2lE + 2 log— , 



(5. 9) 



with 7£ = 0.5772. We note that (5. 8) reduces to the leading thermal effective 
mass uj = M Lead , if the non-leading contributions given through F 2 (u,T) are 
set to zero. The numerical solution of (5. 8) leads to the thermal effective mass 



where p 
M, are 



M = M Lead (l + p), (5. 10) 

0.02223. The massless quark dispersion relations, for the limit k « 



u{k) 



k P 

3 3M v ' 



(5. 11) 



(5. 12) 



where M is given by (5. 10). For the limit k » M, the dispersion relations are 
written as (4. 29) and (4. 30), with M 2 R = M 2 . The inclusion of non-leadine 
contributions leads to increase M^ ea d by 2.223%. We have checked the validity of 
expressions (5. 11) and (5. 12) by comparing the results obtained from them with 
the results obtained by means of the exact numerical solution of the dispersion 
relations. The relative error of the results obtained from (5. 11) and (5. 12) is 
smaller than 0.01% with respect to the exact numerical results. Given the fit 
(5. 10), the solution of (5. 8) can be well approximated by 



UJ 



M. 



Lead 



1 - 2A Z ( lE + log 



M, 



Lead 



2T 



-1/2 



(5. 13) 



Expanding the root square in (5. 13), we can see that the non-leading terms of 
the quark self-energy modify M Lead as oc g ( n) T, with (n = 3, 5, ...), in accordance 



with the conclusions given in [32 



The Lorentz-invariant functions in the Landau gauge are 
a £ (u,k) = ^g 2 [A(0,0) + A 5 (0,0)}, 

b £ (u,k) = ^ 2 [5(0,0) + ^(0,0)], 



(5. 14) 



(5. 15) 



4 The analytical form found for I2, for u> ~ T/2, fits the numerical solutions with a precision 
smaller than 0.01%. 
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where A(0, 0) and 5(0,0) are given by (3. 21) and (3. 22) with M 1 = M 2 = 0. 
The integrals A s (0, 0) and B s (0, 0) are given by (B. 9) and (B. 10) with M 1 = 0. 
Substituting (5. 14) and (5. 15) into (2. 16), we compute the quark dispersion 
relations in the Landau gauge. These dispersion relations coincide with those at 
leading order 0(g s T) in the Feynman gauge, as is shown in Fig. 9. In other words, 
the Lorentz-invariant functions (5. 14) and (5. 15) reduce to the analytical form of 
(4. 21) and (4. 22) with Li p — 4(?g/3. This means that non-leading contributions 
of A(0,0) and 5(0,0) in (5. 14) and (5. 15) are canceled exactly by Ag(0,0) and 
£>q(0,0) contributions, respectively. For k « u, we have analytically proved 
this fact. For the Landau gauge, the gauge dependent parameter p in (5. 10) 
is zero. This result indicates that the one-loop dispersion relations are only 
gauge invariant at leading order 0(g s T) |3|, [33|]. As the thermal imaginary part 
of the one-loop self-energy contains only non-leading terms, it is highly gauge 
dependent. For a consistent calculation of this part, it is necessary to go to the 
imaginary time formalism of the thermal field theory and to use Pisarski's method 
for resumming thermal loops [II, 



We will now study the massive quark case in QCD. The Lorentz-invariant 
functions in the Feynman gauge (3. 18)- (3. 20) contain terms dependent of M q . 
For the case M q « u, it is possible to neglect^ M q in the logarithmic integrals, 
following the same arguments as given in |§. 

We can write (3. 18)-(3. 20) for M q « to as 



a(u>, k) 



2 log 



M n 



ttT 



+ 1E 



(5. 16) 



b(u, k) 



M 2 

bQ(co,k) + A 2 -f 
k 



2 log 



1 



+ 1E 



c(M q ,0) = -AA 2 -I 2 



(5. 17) 
(5. 18) 



where a®(uj, k) and b®(u), k) are the Lorentz-invariant functions for the massless 
case (5. 1) and (5. 2), and I 2 is given by (5. 4). Substituting (5. 16)-(5. 18) in 
(2. 15), taking the limits k << lu, M q « to, and setting k = 0, we obtain the 
following relation: 



1-A 2 F 2 (lu,T) 



~ M 2 Lead 



l-4A 2 F 2 (u,T) 



(5. 19) 



where F 2 (uj,T) is given by (5. 9). This relation for M q = reduces to (5. 8). If 
non-leading contributions manifested by means of F 2 (u,T) are set to zero, the 
relation (5. 19) leads to the leading thermal effective masses M± defined in (4. 13) 
and (4. 14). 

5 We have numerically proved that neglecting the quark mass in the logarithmic integrals, 
for M q < oj/20, leads to relative errors in the integrals smaller than 0.1%. 
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Taking into account the validity of (5. 19) only for M q « M Lead , we obtain 
through a fit of the numerical solutions of (5. 19), that 



M±(M q ) = M Lead (l + p 



2T 



1=F^ 
2T 



(5. 20) 



where p = 0.02223 and a = 0.135. For M q = 0, (5. 20) reduces to (5. 10). The 
dispersion relations given by (4. 13) and (4. 14), with M± defined by (5. 20), con- 
tain all T contributions from the one-loop self-energy. We have also checked the 
validity of these analytical dispersion relations by comparing the results obtained 
from them with those obtained sol- ving numerically the dispersion relations. For 
M q < M Lead /20, the relative error is smaller than 0.01%. 

The Lorentz-invariant functions in the Landau gauge for massive quarks in 
QCD are given by 

a £ (tu,k) = ^g 2 s [A(M q ,0) + A s o (M g ,0)\, (5. 21) 

b £ (co,k) = ^g 2 s [B(M q ,0) + B s (M g ,0)}, (5. 22) 

c £ (u,k) = ^g 2 s C(M q ,0), (5. 23) 

where A(M q ,0), B(M q ,0) and C(M q , 0) are defined by (3. 21)-(3. 23), with 
M 2 = 0. A s o (M q ,0) and B s o (M g ,0) are defined by (B. 9) and (B. 10), respec- 
tively. The quark dispersion relations in the Landau gauge are computed insert- 
ing (5. 21)-(5. 23) into (2. 15). These dispersion relations coincide with those 
obtained at leading order in the Feynman gauge, as is shown in Fig. 10. The 
gauge dependent parameters p and a in (5. 20) are zero in the Landau gauge 
computation; consequently, the massive quark dispersion relations are gauge in- 
variant at leading order in temperature 0(g s T) and in mass 0(M q ). 

The inclusion of the non-leading one-loop contributions changes the leading 
thermal effective masses according to the quark mass. For instance, the change in 
the leading thermal effective masses for the c quark at T = 100 GeV are: 8M + = 
0.0202M+ and 5M_ = 0.0241M_. For the s quark, they are: 5M+ = 0.0220M+ 
and <5M_ = 0.0224M_. Thus, the effect of the sub-leading contributions to the 
dispersion relations is nearly 2% for the light quarks. We have not obtained a 
similar fit of the thermal effective masses as (5. 20) for the heavy quark case. We 
have checked numerically for this case that the sub-leading contributions to the 
thermal effective mass are small in comparison with the leading contributions. 

For the MSM dispersion relations in the unbroken phase, the situation is quite 
similar to the QCD case for massless quarks. For instance, the thermal effective 
masses including leading and non-leading contributions in the Feynman gauge 
are given by (see Eq. (5. 10)): 

Mt R = M u L f ad \l + p), (5. 24) 
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being M^ R ^ ad) the leading contribution to the thermal effective masses in the 
unbroken phase 



M 



u(Lead) 
L,R 




(5. 25) 



with jf,R = f + fin + Jh + fifa-i where the coefficients /„ are given by Eqs. 
(4. 4)-(4. 7). The Lorentz-invariant functions in the Landau gauge for this case 
are presented in the Appendix C. The dispersion relations in the Landau gauge 
coming from substituting (C. 1) and (C. 2) in (2. 13) and (2.14). In a similar 
way as the QCD case, the branches in the Landau gauge are superimposed on 
the Feynman gauge branches at leading order. In consequence, the fermion dis- 
persion relations in the unbroken phase are gauge invariant at leading order. The 
consideration of all terms in the quark self-energy introduces a small gauge de- 
pendence in the dispersion relations. Taking into account (5. 24), the sub-leading 
contribution to the thermal effective mass in the Feynman gauge is 



SMl R = pM u L % ead \ 



(5. 26) 



where p ~ 0.023 for this case. This contribution induces a gauge dependence in 
our numerical computations of nearly 2%. By the mentioned general considera- 
tions this sub-leading contribution is of order ~ g^T in a general gauge. 

The results obtained for massive quarks in QCD can be extended to the MSM 
in the broken case. The quark dispersion relations are gauge invariant at leading 
order in temperature and in quark mass. The non-leading order terms in the 
quark self-energy introduce a small gauge dependence in the dispersion relations. 
We have checked numerically that the thermal effective masses in the broken 
phase M h L R can be written as: 



Mr 



M b + (M q ) L>R + M b _(M q ) LtR 



L L,R 2 " ' 

where taking into account (5. 20), M±(M q ) LtR are given by: 



M b ± (M q ) L , R = M b L % ead) (l + p 



2T 



2 



1-a 



2T 



(5. 27) 



(5. 28) 



with a ~ 0.14, and M b L { R ad) are the thermal effective masses at leading order in 
the broken phase 



M 



b(Lead) 
L,R 



J^2 



f L,R o > 



where the coefficients are given by 



fl,R = / + (!- x z )ffa + (1 - x H )f H + (1 - x w )fZ R . 



(5. 29) 



(5. 30) 



In this equation x^, xh, xw denote the effect of the boson mass contributions 
on the thermal effective masses, being < xz-, xh, xw << 1- This fact permits 
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to understand why the thermal effective masses in the broken phase M b LR are 
smaller than those in the unbroken phase Ml R . The sub-leading contribution to 
the thermal effective masses in the Feynman gauge, taking into account (5. 27), 
is 



b _ iirb(Lead) 



SMI* = PM?X 



2T 



2 



M 

+ (5-31) 



inducing a gauge dependence in our numerical computations about 2% for the 
light quarks. Due the mentioned general grounds, the sub-leading contribution 
to the thermal effective masses is of order ~ g 3 [T — M q ] for quarks, and of order 
~ g 3 [T — M e ] for charged leptons and neutrinos in a general gauge. 

It is easy to prove that the gauge dependence induced in the shift of the ther- 
mal effective masses between the two phases is smaller than the one mentioned 
above. The difference between the gauge dependent contributions to the thermal 
effective masses in a general gauge is of order ~ glM q for quarks and of order 
~ g 3 M e for charged leptons and neutrinos. For the light quark case, the difference 
in the Feynman gauge is: 

5Ml R - 5Ml R ^ - P M b L %^ ^ + (5. 32) 

This gauge dependence is smaller than (5. 26) and (5. 31). 



6 Numerical results 

We now give the numerical results we have obtained for the dispersion relations 
for up quarks (u, c, t), down quarks (d, s, b), charged leptons (e, fi, r) and neutrinos 
(z/ e , v T \ For T = 100 GeV, the following values for the masses in GeV have 
been used: M w = 52, M z = 59, M H = 100, M d = 0.006, M s = 0.09, M b = 2.9, 
M u = 0.003, M c = 1.0, M t = 114.7, M e = 0.0003, M M = 0.069, M T = 1.176 
and massless neutrinos. As we have mentioned, we are specially interested in the 
low-momentum regime where the collective behaviour is present. Adequate boson 
and fermion mass values for each temperature have been used, for the cases where 
the thermal dependence has been studied. The values of the coupling constants 
have been fixed to a s = 0.1, a w = 0.035 for all temperatures. 

In Figs. 3 and 4, we show the branches for quarks s and c respectively in both 
phases. In the unbroken phase, the left-handed branches are above the right- 
handed ones. Although the branches in the broken phase are not associated to a 
specific chirality, the number of branches is also four. While in the unbroken phase 
the left-handed abnormal branch is crossed by the right-handed normal branch, 
in the broken phase there is not crossing between the branches. In the broken 
phase, the branches that now have a hyperbolic form are separated between the 
minimum and the maximum by an energy gap of order ~ Mf. The quasi-particle 
excitations cannot exist in this energy range. 
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The branches for each quark flavour are different in the broken phase and 
clearly the fermionic spectrum described by them is flavour non-degenerate. In 
the unbroken phase, the left-handed branches for the different quark flavours are 
numerically the same in our computation, but a very small difference should exist 
for them. The fermionic spectrum is flavour non-degenerate, but the differences 
for different flavours are very small. The fermionic spectrum described by the 
right-handed branches is also flavour non-degenerate in a similar way to the left- 
handed fermionic spectrum. 

Since temperature and mass corrections are considered, the thermal effective 
masses in the broken phase are smaller than those in the unbroken phase. For 
the s quark case, for instance (see Fig. 3), the thermal effective mass of the left- 
handed branches in the unbroken phase is ~ 51.4 GeV, while that for the 
corresponding branches in the broken phase is ~ 50.0 GeV. The shift of the 
thermal effective masses for the lower branches is also present, but is smaller. In 
our perturbative calculation, this effect can be understood through the presence of 
the massive boson propagator (3. 2) in the fermion self-energy. This discontinuity 
of the thermal effective masses between the two sides of the boundary created 
during the first order electroweak phase transition has recently been noted in 
[ |T3] ] . In this reference, using a numerical simulation of the phase transition in the 
SU(2) Higgs model on the lattice, was found that the thermal effective mass for 
the W and Higgs boson is smaller in the broken phase than in the unbroken one. 

The dependence with T of the top quark thermal effective masses for the 
branches in the broken phase is presented in Fig. 5. Two branches go down 
appreciably as the temperature is lowered (M t increases). For T below 120 GeV, 
these branches can become unphysical as is suggested by this figure.^ The special 
behaviour of the top quark branches is obtained because this quark acquire a large 
mass during the electroweak phase transition. 

The branches for the r lepton are shown in Fig. 6. An energy gap between the 
hyperbolic branches in the broken phase is also present. As in the quark case, the 
thermal effective masses for charged leptons in the broken phase are smaller than 
those in the unbroken one. In the broken phase, the branches for each flavour are 
different and consequently the fermionic spectrum is flavour non-degenerate. In 
the unbroken phase, the branches of a specific chirality are the same for different 
flavours within the precision of the present computation. The fermionic spectrum 
described is flavour non-degenerate, but the differences for different flavours are 
very small. 

The branches for the electron neutrino are shown in Fig. 7. There are only two 
branches because the left-handed chirality of the neutrino. In both phases, the 
branches for different flavours are the same in our computations. However, the 

6 Because of numerical problems, it was not possible to obtain results for the branches for 
T < 155 GeV. 



26 



fermionic spectrum is flavour non-degenerate, in both phases, but the differences 
are very small for different flavours. The thermal effective masses for neutrinos in 
the broken phase are smaller than those in the unbroken one. The T dependence 
of the electron neutrino thermal effective mass is presented in Fig. 8. The thermal 
effective mass is zero for T = 0, as can be concluded from this figure. 



7 Conclusions 

One-loop dispersion relations for the different fermionic sectors of the MSM 
were obtained. The calculation was performed in the real time formalism of the 
thermal field theory, in the Feynman gauge. Numerical results were presented 
for temperatures of the plasma near the critical one of the electroweak phase 
transition. For all the fermions in the broken phase and in the unbroken one, 
the fermionic spectrum described by the branches of the dispersion relations is 
flavour non-degenerate. 

In the unbroken phase, the fermionic quasi-particles can be associated to a 
specific chirality because the chiral symmetry is preserved by the MSM before 
the ESSB. Since the chiral symmetry is violated by the MSM after the ESSB, the 
quasi-particles cannot be associated to a specific chirality in the broken phase for 
the case of massive fermions. The number of branches for quarks and charged 
leptons is four in both phases, while it is only two for neutrinos. 

The thermal effective masses for all fermions in the broken phase are smaller 
than those in the unbroken one. The existence of this shift for the thermal 
effective masses can be understood through the presence of the massive boson 
propagators in the fermion self-energy. A similar shift for the W and Higgs boson 
thermal effective mass has recently been observed in |13|], but through a non- 



perturbative lattice calculation. For quarks and charged leptons in the broken 
phase, an energy gap of order ~ Mf was obtained between the branches that 
have a hyperbolic form. This gap appears as a consequence of the electroweak 
contributions to the dispersion relations. 

The T dependence of the thermal effective masses for the top quark and for the 
electron neutrino was studied numerically. The results for the massive top quark 
suggest that two branches become unphysical in the broken phase for T < 120 
GeV. 

We have studied the gauge dependence of the dispersion relations for the 
massless and massive quark cases in QCD by confronting computations in the 
Feynman and the Landau gauges. We have found for both cases that the disper- 
sion relations in the Landau gauge are equal to those obtained at leading order 
in the Feynman gauge. 

The fermion dispersion relations in the MSM include leading and sub-leading 
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contributions. For the unbroken phase, in a similar way as the massless QCD 
case, the branches in the Landau gauge are superimposed on the Feynman gauge 
branches at leading order. Consequently, the fermion dispersion relations are 
gauge invariant at leading order. The consideration of the non-leading terms 
in the fermion self-energy introduces a small sub-leading gauge dependence in 
the dispersion relations. General considerations indicate that the sub-leading 
contribution to the thermal effective masses, in a general gauge, is of order ~ g 3 T 
for quarks and ~ g 3 T for charged leptons and neutrinos. Particularly, the sub- 
leading contribution in the Feynman gauge is SMI R = pMl { ^ ad) for the quark 

case, where M^^ ad ^ is the thermal effective mass at leading order in the unbroken 
phase (see Eqs. (5. 25)) and p ~ 0.023. The sub-leading contribution induces a 
gauge dependence in our numerical computations of nearly 2%. 

The fermion dispersion relations in the broken phase are gauge invariant at 
leading order in temperature and in fermion mass. The small sub-leading con- 
tribution to the dispersion relation is gauge dependent. On general grounds, 
the sub-leading contribution to the thermal effective masses, in a general gauge, 
is of order ~ g 3 [T — M q ] for quarks, and of order ~ g 3 [T — M e ] for charged 
leptons and neutrinos. The sub-leading contribution in the Feynman gauge is 
5M b LR (M q ) = pM b L (L R ead) [l - where M b L % ead) are the thermal effective 

masses at leading order in the broken phase (see Eq. (5. 29)), and the gauge 
parameters p ~ 0.023 and o ~ 0.14. The gauge dependence induced in our 
numerical computations is about 2% for the light quark case. 

The gauge dependence induced in the shift of the thermal effective masses 
between the two phases is smaller than the one mentioned above. The difference 
between the gauge dependent contributions to the thermal effective masses in a 
general gauge is of order ~ g\M q for quarks and of order ~ g 3 M e for charged 
leptons and neutrinos. This difference for light quarks in the Feynman gauge is: 
5M^ R — bM h L R ~ pM h ^ L R ad '-^7 — cr^r. In our numerical computations is true 
that 5M£ — 5M h L « M'£ — M b L) showing that the gauge dependence in the shift 
is very small. 
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Appendices 



A Calculation of Z-tadpole diagram 

The generic Z-tadpole diagram in Fig. 1 is given by: 

X tad (K) = -i^U E D^(0) J ^Tr^UjSip)}, (A. 1) 

where S(p) is the fermionic finite temperature propagator (2. 1), Z}^(0) is the 
propagator of the Z boson at momentum transfer equal to zero: 

(0) = f^, (A. 2) 

and Ux = (T x — QxS^)L — QxS^R, are the chiral operators of the tadpole 
vertices, with X — E,I. Here, E and / refers to external and internal fermion 
respectively. The expression (A. 1) has the same form as that shown by Notzold 
and Raffel in [2H] for the particular case of neutrino propagation. 



Taking over the expressions of the propagators (2. 1) and (A. 2) in (A. 1) and 
only considering the real part contribution R,eH' tad (K) of the finite temperature 
contribution H' tad (K), we obtain: 



where n(f) is the Fermi-Dirac distribution function given by (2. 6) and Qj = 
(Tf — QiS^)/2. We follow the same procedure as the one used for the evaluation 
of the generic diagrams of Fig. 2. If we multiply (A. 3) by either ^Tor ^, and then 
take the trace of the product, we obtain: 

-Tr(^K ad ) = ^Po-^P)^- W' ~ Mf)n f (p), (A. 4) 

-Tr{iR^' tad ) = J - IP? - Mj)rif(p). (A. 5) 

The argument of the delta-Dirac functions of the expressions (A. 4) and (A. 5) 
have the zeros when p Q = ±(\p\ 2 + Mf) 1 / 2 = e 7 . When we integrate over p 0l we 
obtain the positive and the negative contributions of e/. Immediately, one can 
see that (A. 5) is zero and the integral (A. 4) is seen to be zero after the angular 
integration. It is direct to prove that if we take the trace of (A. 3), the result 
vanishes. As a consequence, the contribution of the Z-tadpole diagram to the 
real part of the fermionic thermal self-energy is zero. This result differs from 
that given in [25|1 , where the tadpole diagram contributes to the real part of the 
neutrino thermal self-energy. This contribution exists because the Fermi-Dirac 
distribution function used there was defined without the absolute value that we 
used in (3. 6). 
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B Real part of the fermion self-energy 
in the Landau gauge 

The thermal propagator in the Landau gauge for massless gauge D £ u (p) and 
massless scalar D £ {p) bosons are: 



PyPu 
9pv + 2 



D £ {p) 



1 



p 2 + ie 



p 2 + ie 



«T 6 (p), 



«r&(p) 



(B. 1) 



(B- 2) 



with Tb{p) defined by (3. 5). The contribution to the self-energy from the generic 
gauge boson one-loop diagram of Fig. 2a is: 



E £ (K) = zg 2 C c J ** D £,(p)>fS<p + lC)>f, 



(27T) 



(B- 3) 



where S(p) is given by (3. 1). Inserting the expressions (3. 1) and (B. 1) into (B. 3) 
and keeping only the real part ReT,' £ (K) of the finite temperature contribution 
J2' £ (K) to the self-energy, it is possible to write: 



Re IT (AT) = ReS'(AT) +Re5E'(K), 



(B- 4) 



where ReE'(iT) is the real part of the finite temperature contribution to the 
self-energy in the Feynman gauge (3. 8) and Re 52' (if) is given by: 



Re6X' = g 2 C c J ^jffcrf + $+ M 



t 6 (p) 



p 2 [(p + K) 2 - Ml 



p 4 



(B.5) 

where the denominators are defined by their principal value. We note that there 
is a principal-value singularity 1/p 2 in the longitudinal boson propagator. Its 
product with S(p 2 ) is defined by —5{p 2 )/p 2 — > S'(p 2 ), where the prime denotes the 
derivative with respect to pi 0. Using this definition and following a procedure 
similar to the preceding case, the traces ^Tr(^Re5T,' ) and ^Tr( r fiRe5T l ' ) are 
calculated. With similar expressions to (3. 15)-(3. 17), the following relations are 
obtained: 



a'(u,k)=g*C e A<(M 1 ,0), 
b 5 (u;,k) = g 2 C c B 5 (M 1 ,0), 



c 5 (uj,k) 



g 2 C c M_ 



2 Mi 



'C(Mi,0), 



(B. 6) 
(B. 7) 

(B. 8) 
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where: 

A s {M u O) = f 
Jo 



dp 1 



47T 2 k 2 



uj 2 - k 2 - M 2 uj 2 - k 2 - M\ 



2p 



+ 



2T 



-e p / T n B (p)[l + 2^L-(p) 



~ (uR (p) + pR + {p)^j - ujp (ujT + (p) + pT (p)) 



'uj 2 -k 2 - Ml {uj 2 -k 2 - M 2 ) 2 



+ 



p_ 

ei 



8A; 

uj 2 + k 2 



16kp 



- + ^e p/T n B {p) 



Lt(p) 



n B {p) 



, , Ltip) - -pS + (p) + oje lP U-(p) + uj 2 pV + {p) 



n F (e 1 ) . (B. 9) 



B 5 o (M u 0) 



dp 1 
o Arc 2 k 2 



uj{uj 2 - k 2 - Ml) 
2p 



~u (ujR'{p)+pR + {p)) - {uj 2 - k 2 )p (ujT+{p)+pT~{p)) 



2\2 



-k 2 - Ml) uj{uj 2 -k 2 - Ml) 



8k 



16kp 



- + ^e p/T n B {p) 
p T 



Ltip) 



n B {p) 



U{U \ k k2) Li( P ) + | P S+(p) + e,{uj 2 " k 2 ) P U-(p)+ 



n F (ei) . (B. 10) 



C(Mi,0) is given by (3. 23) with M 2 = 0. The functions defined in (B. 9) and 
(B. 10) are: 



RHp) 
sHp) - 

T±(p) - 

uHp) - 

V ± (p) - 



{uj 2 -k 2 - Ml 



2\2 



± 



( u * -k 2 -M{ 



2\2 



{uj 2 -k 2 - M 2 + 2ujp) 2 - 4k 2 p 2 {uj 2 -k 2 - Ml - 2ujp) 2 - Ak 2 p 2 ' 

(B. 11) 

{uj 2 - k 2 ){uj 2 -k 2 + Ml) {uj 2 - k 2 ){uj 2 - k 2 + Ml) 



{uj 2 -k 2 + M 2 + 2uje 1 ) 2 - Ak 2 p 2 {uj 2 -k 2 + Ml - 2uje 1 ) 2 - Ak 2 p 2 ' 

(B. 12) 

{uj 2 -k 2 - Ml) {uj 2 -k 2 - Ml) 

{uj 2 -k 2 + M 2 + 2uje 1 ) 2 - Ak 2 p 2 {uj 2 - k 2 + Ml - 2uje 1 ) 2 - 4k 2 p 2 ' 

(B. 13) 

{uj 2 - k 2 ) {uj 2 - k 2 ) 

{uj 2 -k 2 + Ml + 2uje 1 ) 2 - Ak 2 p 2 {uj 2 - k 2 + M 2 - 2uje 1 ) 2 - 4k 2 p 2 ' 

(B. 14) 

{uj 2 -k 2 + Ml) {uj 2 - k 2 + Ml) 



± 



{uj 2 -k 2 + Ml + 2uj £i ) 2 - Ak 2 p 2 {uj 2 -k 2 + Ml - 2uje x ) 2 - Ak 2 p 2 ' 

(B. 15) 

and the functions Lf and Lf given by (3. 13)-(3. 14) with M 2 = 0. 
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The Lorentz-invariant functions in the Landau gauge are given by 

a £ (uj,k) = g 2 C c A £ (M 1 ,0), (B. 16) 

b £ (u,k) = g 2 C c B £ (M 1 ,0), (B. 17) 

c £ (co,k) = ^g 2 C c ^-C(M u 0), (B. 18) 

where 

A £ (M 1 ,0)=A(M 1 ,0) + A*(M 1 ,0), (B. 19) 

B £ (M 1 , 0) = B(Mi, 0) + B S (M U 0), (B. 20) 

with A(M u 0), B(M u 0) and C(M u 0) given by (3. 21)-(3. 23), with M 2 = 0. 
Ag(Mi,0) and ^(M^O) are given by (B. 9) and (B. 10). 



C Lorentz-invariant functions in the Landau gauge 

The chiral projections of the Lorentz-invariant functions for the quarks in the 
unbroken phase are given by: 



a £ {u, k) IF = flA(0, 0) + flA&O, 0), 

R R R 

hi (u, k) IF = flB(0, 0) + flB 5 o (0, 0), 

R R R 



(C. 1) 



(C. 2) 



where the integrals A(0,0) and -8(0,0) are given by (3. 21) and (3. 22), with 
Mi = M 2 = 0. The integrals A s (0,0) and 5^(0,0) are defined by (B. 9) and 
(B. 10), with M 1 = 0. The coefficients are: 



fl = 8 



27ra<- 



t&n 2 9 w 1 2 
27 3 



8 



K+\ 2 K iF , (C. 3) 



11 = 8 



1 + 



tan 2 6> 



w 



27 



/r = 8 



2nas t^olw 



Qjt&n 2 e w + -\] 



SlF, 
8lF, 



11 = 8 



2jras irctw 



Q 2 tan 2 6w 



'IF- 



(C. 4) 

(C. 5) 
(C. 6) 



3 2 

The chiral projections of the Lorentz-invariant functions for charged leptons 



are: 



a{(u, k) = d} L A(0 : 0) + (P L A 5 (0, 0), 

R R R 

h £ L (to, k) = d\B(0, 0) + JlB^O, 0), 

R R R 



(C. 7) 
(C. 8) 
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where the integrals A(0, 0) and -6(0,0) are given by (3. 21) and (3. 22), with 
Mi = M 2 = 0. The integrals Ag(0,0) and B 5 Q (0,0) are defined by (B. 9) and 
(B. 10), with Mi = 0. The coefficients are: 



2 + - + A^ 



2+1 



Ana 



w 



A 2 

tan 2 Ow + -rr 



d R = 4naw 



tan 2 6» 



(C. 9) 

(C. 10) 

(C. 11) 
(C. 12) 

The chiral projections of the Lorentz-invariant functions for neutrinos are: 

af(u, k) = 5^(0, 0) + v 2 A s (0, 0), (C. 13) 

b{ (u, k) = v l B(<d, 0) + ^ 2 ^(0, 0), (C. 14) 

where the integrals A(0, 0) and 5(0,0) are given by (3. 21) and (3.22), with 
Mi = M 2 = 0. The integrals A 5 (0,0) and 5^(0,0) are defined by (B. 9) and 
(B. 10), with Mi = 0. The coefficients are: 



IV 



v = iraw 



2 + i + 



2 + 



(C. 15) 
(C. 16) 
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Figure 1: Generic Z-tadpole diagram 
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Figure 2: Generic gauge and scalar boson diagrams 
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Figure 3: Dispersion relations at T = 100 GeV for the s quark in the broken 
(solid curves) and the unbroken (dotted curves) phases, neglecting quark mixing. 
The upper (lower) dotted curves correspond to left (right) chirality branches. An 
energy gap (~ 0.1 GeV) appears between the minimum and maximum of the 
hyperbolic branches. 
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Figure 4: Dispersion relations at T = 100 GeV for the c quark in the broken (solid 
curves) and the unbroken (dashed curves) phases, neglecting quark mixing. The 
energy gap between the hyperbolic branches is ~ 1.0 GeV. The thermal effective 
masses for the branches in the broken phase are smaller that those in the unbroken 
phase. 
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Figure 6: Dispersion relations at T = 150 GeV for the r lepton in the broken 
(solid curves) and the unbroken (dashed curves). 
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Figure 7: Dispersion relations at T = 150 GeV for the electron neutrino in the 
broken (solid curves) and the unbroken (dashed curves) phases. 
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Figure 8: Variation with temperature of the electron-neutrino thermal effective 
masses for the broken (solid curves) and unbroken (dotted curves) phases 
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Figure 9: Normal branch for the massless b quark in QCD computed for three dif- 
ferent cases at T = 100 GeV. The lower dashed curve is computed in the Feynman 
gauge considering only the leading 0{T 2 ) contribution to the self-energy (leading 
curve). The upper dashed curve corresponds to the results in the Feynman gauge 
considering all the contributions to the self energy (Feynman curve). The solid 
curve is computed in the Landau gauge also considering all the contributions 
(Landau curve). The leading and Landau curves are superimposed. 
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Figure 10: Normal branch for the massive b quark in QCD computed for three dif- 
ferent cases at T = 100 GeV. The lower dashed curve is computed in the Feynman 
gauge considering only the leading 0{T 2 ) contribution to the self-energy (leading 
curve). The upper dashed curve corresponds to the results in the Feynman gauge 
considering all the contributions to the self-energy (Feynman curve). The solid 
curve is computed in the Landau gauge also considering all the contributions 
(Landau curve). The leading and Landau curves are superimposed. 
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